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Abstract. This paper is devoted to the study of the following problem. We have set of 
diffusion processes with absorption on boundaries in some region at initial time t = 0. It is 
required to estimate of number of the unabsorbed processes for the fixed time r > 0. The 
number of initial processes is considered as function of r and r — » oo. 



Consider the N independent diffusion processes which start from different points Xk <E 
Q, k = 1, N in the initial moment of time t = 0. The region Q C R d is open connected 
region and it is limited by smooth surface dQ. All processes are diffusion processes with 
absorption on the boundary dQ. These processes are solutions of the following stochastic 
differential equations in Q with absorption 

d 

dm = a(t, t(t))dt + h(t, m)dwi k) (t) m e R d u) 

bi(t, x), a(t, x) : R+ x R d -> R d . 
with an initial condition: £(0) = 7& €E Q. 

Here the W {k) {t) = {wf\t), l<i<d), 1 < k < N are independent in totality d- 
dimcnsional Wiener processes. 

Thus, these processes have the identical diffusion matrices and shift vectors , but they 
have different initial states. 

We will be interested by distribution of number of processes yet not absorbed by 
boundary dQ in the moment of time r. 

This task was offered in [1] as the mathematical model of cleaning of gas from particles 
(dust, microbes and ect.). Cleaning consists in pass of gas with speed v through the pipe 
of length I and with permanent section Q U dQ . Walls of pipe absorb the particles. 
The estimation of efficiency of cleaning can be reduced to the solution of this problem 
at t = I /v. In supposition, that particles evolve independently from each other and their 
movement can be presented as the solutions of (1). 

We make the following assumptions: a number and positions of particles are defined 
by a given determinate limited measure N(B, r) at the initial time. Thus N(B, r) is 
equal to number of points xu in the set B and N = N(Q, t) < oo for fixed r > 0. 
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Note, the case when N(-, r) = N(-) there is random Poisson measure on a circle 
C = {(x, y) : x 2 + y 2 < r 2 } was considered in [1]. 

There is exact formula of distribution function for number of remaining processes in 
this article. However, this formula consists diffucult computed functions. The present 
article is devoted to case when initial number of processes N(Q, r) depends on the final 
time t: N(Q,t) — > oo when r — > oo. We shall obtain conditions of such dependence 
which leads to simple limit distribution function of number not absorbed processes. 

Suppose that the region Q is bounded and boundary dQ is Lyapunov surface C^ 1 '^ 
[2] . We will consider the following case 

a(t,x) = (0, . . . ,0), bi(t,x) =bi = (bii,...,bid), 1 < i < d. 

d 

We will define matrix a = B T B, B = (hj), 1 < i,j < d a = (cry ), 1 < i, j < d and 
differential operator A : a ij dx dx ■ ■ CT be a matrix with the following property 

l<i,j<d * 3 

^ <7i ZiZj > H\z\ 2 . 

l<i,j<d 

Here /i, there is fixed positive number, and z=(zi,--- ,Zd) there is an arbitrary real 
vector. 

This operator acts in the following space 

H A = {u : u e L 2 {Q) n Au e L 2 {Q) n u{dQ) = 0} 

with inner product (u, v)a = {Au, v). Here (,) is inner product in L2(Q). The 
operator A is positive operator [2]. It is known that the following eigenvalues problem 

Au = -Au, u(dQ) = 
has infinity set of real eigenvalues \ — > oo and 

< Ai < A 2 < • • • < A s < • • • . 
The corresponding eigenfunctions 

/ll 7 ■ • • 7 /l«l 7 ' ' ' 7 JSI 7 ' ' • 7 fsU s , ' ' ' 

form complete system of functions both in Ha and L 2 {Q) := {u : u e L 2 (Q)r\u(dQ) = 0}. 
Here number nk is equal to multiplicity of eigenvalue A^ . 

We will denote by t](t) the number of remaining processes in the region D at the 
moment r. 

We will also assume that a-additive measure v is given on the £„- algebra sets from 
Q, v{Q) < oo. All eigenfunctions : Q — > i? 1 arc (E, y ,Sy) measurable. Here Sy is 
system of Borel sets from i? 1 . Let ^> denote the weak convergence of random values or 
measures. 

Let us denote by v T {-) the measure 

v T {B)=ex V {- T -\ 1 )N{B,T). 

where BeE„. 

By definition of measure v T { ), we have 

j f exp(-§Ai), if x = x k , k = 1,- • • ,N(Q,t) 

I 0, otherwise. 
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Theorem 1. Suppose N(-,t) satisfies the condition 

M-) => "(■)■ 

r — >oo 

Then w(t) => n if t — > oo where n has Poisson distribution junction with parameter 

n 1 

a = J F(x)dv(x) and F(x) = J2 fu(x)c u , c u = J f u {x)dx. 
Q »=i Q 

Proof. Consider the following initial-boundary problem 

du 1 x - d 2 u 

M = 2^° ij dx~dx~ j > XeQ; 

l<zj<rf J 

u(0, x) = 1, x e Q 



u(t,x) = if xe3Q, t>0 (2) 

It is known [3, sec. VIII], that u(t,x) is equal to probability to remain in the re- 
gion Q at the moment t for a diffusion process which was in a point x at the initial 
moment ( £(0) = x £ Q). 

Introduce indicators 

f 1, if k — th particle belongs to Q at the moment r 

X/s(t) = < 

[0, if fc — th particle absorbed by the moment r. 

These indicators are mutually independence by assumption. Thus the following rela- 
tions are correct 

N(Q,r) 
fe=l 

N(Q,t) N(Q,t) 

Es"W= [] Es^= [] («*(r,7fc) + (l-«(r,7fc)); 
fe=i fe=i 

AT(Q,t) 

lnSs"( T '= ^ In(l-u(r,7fc)(l-a)). 

Here < s < 1. As < u(t, -fk) < 1 then we have the following inequality from the 
last 

N(C,r) N(C,t) 

\]nEs^- «(r, 7 fc)(l-a)l<« £ U V>7*), (3) 

fe=i fc=i 

Here a < oo. Define the value of u(r, 7fe). 

We shall assume that system of functions {fy (x),i> 1,1 < j < Ui} is orthonormalized 
with respect to space L^Q). The ordinary argumentaion (see, for example, [2, sec. 22]) 
loads to the definition of solution of problem (2) of the form 
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i — i n k 

u{t,x) = ^exp(--A fe )^c fcj / fej (a;) 
fe=i j=i 

where coefficients are equal to coefficients of decomposition of initial value (unit) 
by system of functions fay. = J fij(x)dx. The Parseval - Steklov equality is true for 

Q 

these coefficients: 

oo n k 

EE c « = m- (4) 
fc=i j=i 



Further 



N(Q,t) 



E u ( T '7fc)= / u(T,x)dv T {x) = 
Q 

J F(x)du T (x) + j s T (x)dv T (x), (5) 



Q Q 



here 



s T (x) := E CX P ("^( A fe _ A i)) E Cfe ^'( x )- 

fe>2 j=l 

As function F(x) is continuous and bounded function on the Q then under the condi- 
tion of the theorem, we obtain [4] 

J F(x)du T (x) — ► J F{x)du{x). 
Q Q 

In order to estimate of s T (x), we will give the result from monorgaphy [5,Thm. 17.5.3] 
Consider the sums of eigenfunctions of the form 

_ n k 

e(x,\)= E E/«(») 

Xk<Xj=l 

then 

sup \fe(x, A) < CA 5 . 

The asymptotic of eigenvalues Afc under fc — > oo is defined by the following eniqualities 
[2, sec. 18] 

Cifc 3 < Afe < c 2 fc d , where ci, c 2 = const. 

The last and (4) and Caushy-Bunyakovskii inequality lead to the following convergence 
under r — > oo 
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\S T {X)\ <E eX p(~^( A fe" A l)) \ E' 
k>2 \ j=l 



'kj> 



<C^a| exp(-^(A fe -A)) 



k>2 



<C\J24eM-T(\k-\i))\ (EE 

ifc>2 j = l 



c kj 



0. 



(6) 



, k>2 



Combining the condition of the theorem and (6), we conclude that second summand 
of right part of (5) converges to zero when r — > oo. 

Now we shall estimate of the right part of inequality (3). 
By definition of measure v T we have 

N(Q,r) 

R T := ^2 u2 ( T >7fc) = / [Ji(r,x) + J 2 {t,x) + J 3 (t,x)} dv T {x). 

I i J 



fc=l 



Here 



E ci ^y( a 

J 2 (r,a;) := 2 ^ cy/i^x) E exp(--A fe ) E c ^ /iW fa); 



k>2 



= 1 



J 3 {t,x) := exp(--X 1 )sl(x). 
Observe that we have F(x) > 0, x £ Q, because the function u(t, x) > for all 

t>0, xeQ. 

Let M = ma,xF(x). Applying (6), we get under r — > oo 

a;€Q 

|Ji(r,a;)| < exp(-^Ai)M 2 - 0; 

|J 2 (r,a;)| < 2Mexp(-^Ai) sup |« T (a;)| -» 0; 

|^3(r,a;)| < exp(-^Ai) supsj(i) -> 0. 
^ xeQ 

These inequalities and theorem's condition guarantee the convergence of R T to zero 
under r — > oo. 

The proof of theorem is complete. 



Example. Now we shall investigate particular case of the general problem. Here we 
can calculate relevant values of normalizing function and on the other hand this case 
may be represent the first approximation of real situation. 
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Consider circle domain Q in E 2 : x 2 +y 2 < r 2 . Assume that the diffusion particles start 
from point (xk,yk) G Q at the moment t = 0. The movement of particles is described by 
the following stochastic differential equations 

2 

d£(t) =Y,bidwi(t) (7) 

i 

£(0) =£o = (a*, y fc ). 

where &i = (cr, 0),&2 = (0, a) and W(t) = (wi(t),i = 1,2) be a 2-dimensional Wiener 
process. 

Assume that the equation (7) defines a diffusion process with absorption on the bound- 
ary dQ = {(x, y) : x 2 + y 2 = r 2 }. 

In what follows, the Jo(x), J\{x) are Bessel functions zero and first order. It are 
defined as the solutions of next equations 

d 2 y 1 dy , n 2 . 

— - H + (1 ) = 0, 

dx 2 x dx x 2 

y{x ) = 0, (x = VAr); |y(0)| < oo; 
under n = and n = 1. 

The value of is equal to to- th root of equation Jo{p.) = [6,7]. 
We will use the symbol rn.es (•) to denote the Lebesgue measure. 
Put 



/(r) := cxp 



2 \ r 




Suppose N(-,t) satisfies the condition 

lim N(B, r)/(r) = mes(B), B e S mes (C). 

The initial position of fc-th particle in this case has form 7^ = (xk,yk) 
The value of u(r, 7^) is defined as value of u(t,x,y) in this point. Here u(t,x,y) be a 
solution of the following initial-boundary problem 

du a 2 ( d 2 d 2 \ , s „ 



dt 2 \9y 2 9a; 2 / 

u{0,x,y) = l (x,y)eQ; 

u(t,x,y) = if (x,y)edD, t>0 (8) 

According to results from book [6, sec. VI] the solution of (8) in the point (r,7fe) can 
be expressed in form 



9 t I ..(0) v^f+g 
-•Ai I /' 



rf 1 
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Let us compute the parameter a. In this case it is convenient to decompose of circle 

C by concentric circles for constuction of integal sums of integral J F{x)v(dx) [8]. 

C 

Define this partition 

K ni = l(x,y) eC: - < ^fx T +tf< ?ll±l21 0<*<n-l. 
I n n J 

Now mes(K ni ) = g(^) - g(±) where g(p) = nr 2 p 2 , < p < 1. 
Finally, the parameter of Poisson distribution is equal to 

l 

a = 2 (^ 0) Ji(/4 0) )) ^ 2^r 2 J J {) {^ ] p)pdp = n 

o 

We used the following known relation aJo(a) = [aJi(a)]' [6, sec. VI] for calculation of 
the last integral. 
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